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Abstract 

< i 

An approximate method is proposed for investigating complex- 
^ . temperature properties of real-dimensional spin-glass models. The 

method uses the complex-temperature data of the ferromagnetic model 
C on the same lattice. The universality line in the complex-temperature 

space is obtained. 

g 
o 
,o, : 1 Introduction. 

Investigating of physical quantities in the complex plane can reveal new and 
(T) ■ unexpected effects. More than four decades have passed since Dyson [|TJ 

considered imaginary charges in electrodynamics. Having a = e 2 /he < 0, 



he got that this theory should be essentially unstable, and expansions in a 



are, at best, asymptotic. The consideration of complex magnetic field by Lee 



and Yang and complex temperatures by Fisher opened new and effective 
method for investigating phase transitions. Later the method found large 
variety of applications in obtaining the order and type of possible phase 
transitions [ 15| , critical indices ||[|13], and recently has been connected with 



experiments on magnetization fT6|| . 



The method can be even more useful when investigating phase transitions 



in disordered systems, since many well-developed analytical and numerical 



methods are not suitable here. This point of view was first provided in [] I 



where so-called Graffiths' singularities was discovered in statistical systems 
with random-fields. Recently, a strongly-frustrated but nonrandom system 



was investigated in complex temperature plane [0 , and a continuous spectra 
of phase transitions was obtained. 



The consideration of spin-glasses in complex-field and/or complex-temperature 
plane was started in ||||. Namely, Random Energy Model (REM) was in- 
vestigated, which is the simplest but typical representative of spin-glasses 
||]. This has been continued in where the more physical dilute (finite- 



connectivity) REM |T(J [11] was investigated. In particular, we have pointed 
out also how the Dilute Generallized Random Energy Model (DGREM) 
should be considered in the complex planes @. This seems very impor- 
tant, because DGREM is nothing else, but the quite accurate approximation 
to the real-dimensional (non-mean-field) Edwards- Anderson (EA) spin-glass 
model. We believe that the results obtained with DGREM can be relevant 
for that model, and will help to clarify the properties of spin-glasses in finite 
dimensions. Recently, the lower critical dimension of the EA model has been 
estimated in this fashion || . This line of research will be continued in the 
present paper. We shall give the density of the partition function's zeroes, 
and discuss in details its applications to the real-dimensional EA models. 

2 Dilute REM 

We shall start with repeating some facts about diluted REM and GREM. 
The dilute p-spin glass model is described by the following hamiltonian 

aN 

H = - Y. 3ii..i p s il ..s ip , (1) 

(n<i2--<jj>)=l 

where only aN couplings ji x ..i p are non-zero 1 > i > N, and Sj = ±1 are 
Ising spins. At high temperatures the system is in the paramagnetic phase 
with free energy 

\nZ = aN\ncosh((3) + N\n2 (2) 

At critical temperature T c = 1//3 C a phase transition occurs into the SG phase 

ag(ta.nh/3 c ) = 1 (3) 

where the function g(x) is defined as 

g(x) = l -{l + x) ln(l + x) + ^(1 - x) ln(l - x)\ (4) 



Below this temperature we have 

\nZ/N = aN0y(\n2/a) (5) 

where y is inverse function to g(x) 

\{l + y) hi(l + y) + ^(l-y) ln(l -y)]=x (6) 

3 Dilute GREM 

Let us consider the diluted version of GREM, with infinite levels of the hier- 
archy. Now we have an infinite chain of DREM-s. At the interval [v, v + dv] 
we have Ns'(v)/ In 2 spins with Nzdv couplings. Function s(v) is monotonic, 
s(0) = 0,s(l) = In 2. Similarly to the case of dilute GREM at real T we 
found: 

0F r^W s'(vn) 

-^- = z(l-v c (0))lncosh0 + (ln2-s(v c ) + z0 dv y(^^) (7) 

JM Jo z 

where v c (0) is defined from the equation 

tanh(/3) = y{s'{v 2 )/z) (8) 

For the case of Edwards-Anderson model placed on d-dimensional hyper- 
cubic lattice 

u , s i n S(-vdN) 

z = d ' v = ~Wd> s ^) = ln2 - N w 

here U is energy, and S(U) is entropy as function of the energy for corre- 
sponding ferromagnetic Ising model. It is easy to derive from the definition 
of temperature 

So there is a connection between and v. For a given we find energy of 
corresponding ferromagnetic model, and then calculate 

v = -U(~0(v))/(Nd) (11) 

We obtain for the free energy 

0F f"c(J3) 

-'^- = z(l-v c (0))\ncosh0(\n2-n(v c (0)))+z0 dv y((3(v )) (12) 

JM Jo 



4 Complex temperatures 



Let us consider now the case of complex temperatures. Now we have 3 phases 
for REM. 



PM: 



SG: 



LYF: 



—— = aRe In cosh(/3) + In 2 (13) 

]UZ alhv(l/a) (14) 



N 



\nZ a „ . , ,_, . . In 2 . . 

— = -Re In cosh{2^) + — (15) 

There are boundaries between SG- PM, SG - LYF ,PM - LYF. 
For the PM-LYF we have a line 

. 2 V"-l-tanh 2 A 

Sm (A) = 2V a( l-tanh^ 1 ) (16) 

This line begins at /3i — and is contined till intersection with SG-LYF line. 
For it we have 

A = A), oo > p 2 > p 2c (17) 

where fac is defined from the intersection with another line, and j3\ c from the 
equation 

a/2 In cosh(2/3 lc ) + In 2 = afi lc y{\n 2/ a) (18) 

Then the third line PM-SG is defined from the equation 

aite In cosh /? + In 2 = /3 lc ay(m2/a) (19) 

Let us vary the parameter a. We can construct some universal line for the 
critical {3 c ,i,Pc,2- If we define function 

fU A _ Ml + t) 2 -ln(l-t) 2 [l-(l-t 2 ) 3 ] 

/lS ' tJ " ln(l + t)/(l-t) W 

then for the critical t = tanh(/5i), s = sin 2 (/3 2 ) we have an equation 

i(l+/( S , t)) ln(l+/( S , t))+I(l-/( s , t)) ln(l-/( S , *))] = ln(l+t 2 )/[l-(l-t 2 ) S ] 

(21) 



5 DGREM at complex T 



Now we have 

0F 

- ^— = d(l - v 2 )Re In cosh/? + (In 2 - s(v 2 )) + d(v 2 (0) - v 1 {0))/2 In cosh 20 x 

(22) 

+[s{v 2 ) - s{v 1 )}/{2) + d0! H dv o y{0{v o )) 

Jo 

The values of v%, v 2 are defined from the extremum condition. Integrating by 
parts in the last term we get 

0F 



N 



d(l-v 2 (0))Re\ncosh0+(\n2-s(v 2 ))+d(v 2 -v 1 )/2\ncosh20 1 (23) 
+(s(v 2 (J3) - s(«i(/5)))/2 - dPx f 1 d^o^r^ + dv&yfa) 

1-2/ 

where y as a function of 0q is defined from the equation 

V = 9^) (24) 

function Vq{0) is defined from the equation 

v o 0) = -U0)/(Nd) (25) 

Here U(0i) is the energy of ferromagnetic model on the same lattice. The 

variables V\,v 2 we can define from the saddle point condition of free energy. 

Let us collect all these results together. We consider some point of hierarchy 

v. At inverse temperatures 0\ it could be in SG or PM phase. If it was in 

SG phase, then it stays there while adding imaginary 2 . 

When at real model was in PM phase there are three possible scenarios of 

behavior: 1) stay in PM: 2) become SG phase; 3)become LYF. We can look 

under another point. What happens with our system, when we vary a at 

fixed values of p\, 2 1 

If our point at complex space is under the line |2I|, then our system could 

be in either SG or PM phase. Above that line all 3 phases are allowed for a 

different parts of hierarchy. 

This line stays for any version of dilute GREM models. Thus it means some 

universality. It will be interesting to find similiar universality classes in other 

mean field models. 



6 LYF zeros at the border of PM SG phases 

Let as first consider the case, when our system is under the line (21) and 
there is not LYF phase. We have an expression 

In Z/N = d(l - v c )Re In cosh (3 + In 2 - s(v c ) + dfc f' dv y((3(v )) (26) 

Jo 

where the value of v c is defined from the saddle point condition 

dRe\ncosh(3 + s'(v c ) = d/3 iy (fiM) (27) 

For the laplacian we have an expression 

dV (f£ d 2 f dv c 2 ( _duc_ ,. d 2 f dv c d 2 f dv c 

K d(3l d(3l } d 2 v c [V d^ } [ d[3 2 ) 2l) dvjdp^dfa' dv c dp 2 V df3 2 } l ' 

To calculate density we need in expressions 

it \ A o i mo nt \ d2 t5\ ^elncosh/5 + /3i 
s (v) =d(3, y = 1/(3, s (v) = — -~-, y{l3 x ) = (29) 

c((3) Pi 

here c is specific heat of ferromagnetic phase. Eventually: 

7Tp(A,ft)= (30) 

d 2 tanh/5isin 2 /5 2 n2 

— {[(y-tanh^i 



s"(-l + /3ij/') cosh/^sin^ + cos^tanh 2 /^ 
r sin/3 2 cos/J 2 2 



7 Three phases 

When our point in complex (f3±, f3 2 ) space is over the universality line, we have 
the follwoing expression for the laplacian: 

(^2+^2) ^ f^i/3i + ^fpiv 1 V l,l3i + fv 1 v 1 ( v l,p 1 ) + 2 f/3iv 2 V 2,/3 1 +fv 2 v 2 ( V 2,/3 1 ) (31) 



For the i>i we have an equation 

fM + rf lncosh2ft=ft!/( £M )d (32) 

Having values i>i,i>2, wwe can define the values of corresponding /3 1? /? 2 by 
means of formula (25). Equation for the v<i- 

s'(t;2)+rflncosh/3i+rf/21n(cos 2 /32+tanh 2 /5isin 2 /32) = s'(u 2 )/2+d/21iicosh2/3i 

(33) 
Eventually we have for zeros density: 

2z(,j 2 - vx) <f(-tanh2ft+j/(fl)) 2 

"^' M = ^^m — cw - 1/2)^) (34) 



2 or (sin /9 2 cos # 



2j i2 



s" (v 2 ) [(cosh /3f (cos /31 + tanh/?! sin # 



, r+ ha + to sm ^smh/^i 2 

+ tanh2pi - tanhpx — — r- -= ■ — — - — : — -j- } 

L cosh/3 1 3 (cos/5| + tanh/5ism/3|) J J 

One can use formulas (30) and (34) to calculate density of zeros, having 
the data of the ferromagnetic model at the same lattice. 



8 Conclusion. 

This paper devoted to the approximate method, which allows to investigate 
the zeroes of the statistical sum for Edwards- Anderson models in real phys- 
ical dimensions. The key point of the method is in the using the rich phase 
spectrum of Dilute Generalized Random Energy Model. We obtained a uni- 
versality line in the phase diagram of the model. It hardly controls the cor- 
responding complex-temperature behavior. The various notions that entered 
our discussion give good hope for the applicability of Random Energy-type 
models to realistic systems. 

Acknowledgements 

D.B. Saakian is grateful to Fundacion Andes grant c-13413/1 for financial 
support. A. Crisanti, S. Kobe and W. Janke are acknowledged for kind 
hospitality. 

E.E. Vogel is grateful to Fondecyt 1990878. 

7 



References 

[1] F.J. Dyson, Phys. Rev. 85 (1952) 631 

[2] C.N. Yang and T.D. Lee, Phys. Rev. 87 (1952) 404; Phys. Rev. 87 
(1952) 410. 

[3] M.E. Fisher, in Lectures in Theoretical Physics, ed. by W.E. Brittin, 
(University of Colorado Press, Boulder, 1965), vol. 7c. 

[4] B. Derrida, Phys. Rev. Lett. 45 (1980) 79. 

[5] B. Derrida, Physica A177 (1991) 31-37 

[6] C. Moukarzel, N. Parga, Physica A177 (1991) 24-30 

[7] D.B. Saakian, |cond-mat/ 9709238; ; JETP Lett. 67 (1998) 5 

[8] D.B. Saakian, |cond-mat/9905123| . 

[9] B. Derrida, M.R. Evans, E.R. Speer, J. Stat. Phys., 156 (1993) 

[10] D. Dominicis, P. Mottishow, J. Phys., A20 (1987)L267 

[11] A.E. Allakhverdyan, D.B. Saakian, Nucl. Phys V.B498 (1997)604. 

[12] A.E. Allahverdyan, A.S. Ananikian, S.K. Dallakyan, Phys. Rev. E, 57 
(1998) 2452. 

[13] C. Itzykson, R.B. Pearson, and J.B. Zuber, Nucl. Phys. B, 220 (1983) 
415. 

[14] R.B. Griffiths, Phys. Rev. Lett., 23 (1969) 17. 

[15] K.-C. Lee, Phys. Rev. Lett. 73 (1994) 2801. 

[16] Ch. Binek, Phys. Rev. Lett, (in press) 



